Acoustic horizons for axially and spherically symmetric fluid flow 
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We investigate the formation of acoustic horizons for an inviscid fluid moving in a pipe in the case 
of stationary and axi-symmetric flow. We show that, differently from what is generally believed, the 
acoustic horizon forms in correspondence of either a local minimum or maximum of the flux tube 
cross-section. Similarly, the external potential is required to have either a maximum or a minimum 
at the horizon, so that the external force has to vanish there. Choosing a power-law equation of 
state for the fluid, P oc p", we solve the equations of the fluid dynamics and show that the two 
possibilities are realized respectively for n > — 1 and n < — 1. These results are extended also to the 
case of spherically symmetric flow. 



Sonic horizons can be generated by the motion of a classical fluid l]. The simplest case is represented by a fluid 
with space-dependent velocity flowing trough a tube with variable section (Laval nozzle) 3 • If ttic fluid flows in the 
direction where the tube becomes narrower, the fluid velocity v will increase downstream. Eventually, a point will 
be reached where v is equal to the local speed of sound c and v > c beyond this point. An acoustic perturbation 
generated in the supersonic region cannot reach the upstream subsonic region. We have generated a sonic horizon, 
the acoustic analogue of a spacetime horizon. The downstream subsonic region behaves, as long as only acoustic 
disturbances are concerned, as the acoustic analogue of a black hole, a dumb hole. 

Starting from the original proposal of Unruh , acoustic (and other condensed matter analogues) black holes have 
been used to describe various kinematical aspects of general relativity such as event and cosm olog ical horizons, field 
theory in curved spacetime, Hawking radiation etc. |1 S S E II i, S ES IHIilllllllllllllllillllli- Recently, 
the fluid/gravity analogy as been extended also at a dynamical level. It has been used to gain information about 
Einstein's equations governing the dynamics of a spherically symmetric gravitational black hole [20I l2l | . The lack of 
direct experimental tests has always been an obstacle for the research on black holes. Investigations on "artificial" 
black holes could represent a useful bridge to make black holes experimentally more accessible. 

Although conceptually very simple, the realization of acoustic horizons requires fine-tuning of the external potential 
acting on the fluid, of the tube profile and of the initial fluid velocity [a,l2ll22l- the case of the Laval nozzle, it is 
widely believed that the horizon must form exactly at the narrowest part of the nozzle However, it is 

not easy to find necessary and sufficient conditions for the formation of the horizon. In the general case, the equations 
governing the dynamics of the fluid are very difficult to solve. 

In this paper we investigate the formation of acoustic horizons in the case of a inviscid fluid with stationary, axially 
or spherically symmetric flow. When considering axi-symmetric fluid motion we will discuss separately the case of 
a flux tube of varying section and the presence of external forces. We show that, differently from what is generally 
believed, the acoustic horizon forms in correspondence of either a local minimum or maximum of the flux tube cross- 
section. Similarly, the external potential is required to have either a maximum or a minimum at the horizon, so 
that the external force has to vanish there. Choosing a power-law equation of state for the fluid, P oc p", we solve 
the equations of the dynamics and show that the two possibilities are realized respectively for n > — 1 and n < — 1. 
Moreover, we will determine the range in the initial conditions for which the horizon effectively forms. 

Axi-symmetric flow 

Let us first consider axially symmetric, stationary, inviscid fluid motion, which is constrained on a pipe. Indicating 
with X the coordinate along the symmetry axis, the fluid is described by its velocity v(x), density p{x) and pressure 
P{x). The flux tube is characterized by its cross-section A{x) and we will assume that external forces, characterized 
by the potential ^(x)), act on the fluid. Because of the symmetries of the problem all the parameters are functions 
of the coordinate x only. The equations governing the dynamics of the fluid are the continuity and Euler equations 
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We will restrict ourself to the case of isentropic fluid flow P = P{p), which gives for the local speed of sound c = 
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FIG. 1: Flow situation (on the left) and qualitative behavior of the parameters (on the right) for a converging- diverging nozzle. 
The black, grey and dashed lines represent the tube cross-section A, the fluid velocity v and the speed of sound c, respectively. 
The vertical axis indicates the position of the horizon, c and v are normalized to their horizon values. 



Wc will discuss separately the case of a flux tube of variable cross-section and homogeneous external potential and 
that of a flux tube of constant section and variable external potential. 

Flux tube with variable cross-section 

Using Eqs. |Q with -0 = const., one easily obtain the well-known nozzle equations |RE^ 
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where M = v/c is the Mach number. The nozzle equation ^ can be used to discuss the conditions for the formation 
of sonic horizons, i.e of a surface with M = 1 separating a subsonic (M < 1) from a supersonic (M > 1) region. 
Assuming that v,v'{x) = dv/dx,A are everywhere finite and nonvanishing, one can easily see from Eq. ^ that, 
the if a horizon forms, this necessarily happens in correspondence of a local extremum of the tube cross-section A. 
Conversely, if the cross-section A{x) has a local extremum at a point Xh necessarily a sonic horizon must form at this 
point. Depending on the sign of v' (x) we have two different situations. 

(i) The fluid velocity grows monotonically passing from the subsonic to the supersonic region. In this case the 
cross-section A has a local minimum at the horizon. This corresponds to the usual situation. The flux tube has the 
shape of the Laval nozzle, a converging pipe where the fluid is accelerated, followed by a throat and by a diverging pipe 
where the fluid continues to accelerate. The corresponding flow situation and qualitative behavior of the parameters 
are shown in Fig. (Q. 

(ii) The fluid velocity decreases monotonically passing from the subsonic to the supersonic region. In this case the 
cross-section A has a local maximum at the sonic horizon. This flow situation is very unusual and until now has not 
be considered in the literature. The flux tube has the form of a diverging pipe where the fluid is decelerated followed 




FIG. 2: Flow situation (on the left) and qualitative behavior of the parameters (on the right) for a diverging-converging nozzle. 
The black, grey and dashed lines represent the tube cross-section A, the fluid velocity v and the speed of sound c, respectively. 
The vertical axis indicates the position of the horizon, c and v are normalized to their horizon values. 
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by a mouth and by a converging pipe where the fluid continues the deceleration. The corresponding flow situation 
and quahtative behavior of the parameters are shown in Fig. 

Notice that if the tube cross-section is given as a function of the fluid velocity, we can write the nozzle equation in 
the equivalent form 

^ = -(A^^-l). (3) 
av V 

This equation relates the appearance of an horizon directly with the existence of local extrema of A{v) without the 
need of further conditions on the behavior of the space derivatives of v. We will make use this equation when discussing 
the formation of horizons in a fluid with spherically symmetric flow. 

Until now our discussion has been focused on the nozzle equation (0) . The previous results are the whole information 
we can extract from this equation. There are two main points on which Eqs. Q give no information. First, they do 
not say to us which of the two possibilities (i), (ii) is actually realized. Second, our results are based on the rather 
strong assumption v'{x) 7^ oo, 0. In general its validity depends on both the dynamics and the initial conditions 1^. 
In order to answer to the previous questions one has to choose an equation of state for the fluid and to solve Eqs. 
To be more concrete we will consider a fluid with a power-law equation of state, 

P=-p-. (4) 

n 

where a,n ^ are real constants. This equation of state describes almost all physically interesting fluids: perfect 
fluid (n = 1), Bose-Einstein condensate (n = 2), Chaplygin gas {n = —1). Even in the simple one-dimensional case 
under consideration explicit solution of Eqs. CD,Q cannot be found. However, one can flnd solutions in implicit 
form, writing p, c, A as a function of the fluid velocity v. The solution of Eqs. for a fluid with equation of state 
given by Eq. Q and n ^ 1 reads. 



pW=(^(«--^))"", c\v) = ^(a~v% A{v) = l-, (5) 

where a, (3 are integration constants determined by the initial conditions. Formation of the sonic horizon requires 
c{vh) — Vh, which yields for n 7^ —1 

V n + 1 

where vu is the fluid velocity at the horizon. Formation of the sonic horizon is possible in the following range of 
variation of the parameters. For n > 1, a > and v < yja = \J(n — l)/(n -I- 1) vn- For — 1 < rt < 1, a < 0. For 
n < —1, a > 0, V > ^/a. The Chaplygin gas {n — —1) represents a limiting case. Formation of the horizon requires 
a = 0, which in turns implies A — const and v = c, identically. 

Using Eqs. © we can also determine which of the two possibilities (i), (ii) is actually realized. Differentiating 
two times A{v) we get {dA/dv)]^^ = and > (< 0) for n > — 1 (n < —1). Thus, the horizon 

forms in correspondence of a local minimum (maximum) of the function A(v) for n > —1 (for n < —1). Moreover 
d'^A/dv^ is everywhere nonzero excluding the possibility of a flex point. To infer about the behavior of A{x), once 
the behavior of A{v) is known, we need just to use the trivial identities dA/dx — {dA/dv){dv/dx), d^A/dx"^ — 
{d'^A/dv'^){dv/dx)'^ + {dA/dv){d'^v/dx^). When dv/dx is finite and nonvanishing, local maxima (minima) of A{v) 
correspond to local maxima (minima) of A[x). We have therefore shown that in the case of a fluid with a power-law 
equation of state |0J, the acoustic horizon forms in correspondence of a minimum (maximum) of the cross-section 
when n > —1 {n < —1). The Chaplygin gas {n = —1) represents a limiting case for which the cross-section must be 
constant and v — c everywhere along the tube. 

Flux tube with constant cross-section and non-homogeneous external potential 

This case can treated similarly to the previous case. Setting A — const in Eqs. (^3) one derives an equation similar 
to the nozzle equation in which A is traded for 7p 



dx dx^ ' 



(7) 
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FIG. 3: Qualitative behavior of the parameters for a flux tube of constant section and varying external potential. The figure 
on the left (right) corresponds to case I, (II), respectively. The black, grey and dashed lines represent the external potential tp, 
the fluid velocity v and the speed of sound c, respectively. The vertical axis indicates the position of the horizon, c and v are 
normalized to their horizon values. 



The main difference between Eq. ^ and is the change of sign in the right hand side. This means that, differently 
to what happens for v'{x) and A'{x) in the nozzle Eq. (O, v'{x) and i/j'ix) have the same (opposite) sign in the 
subsonic (supersonic) region. 

Again, assuming that v and v'{x) are everywhere finite and nonvanishing, one can easily show using Eq. Q that 
the horizon must form in correspondence of local extrema of the external potential "0. This gives a a null (external) 
force condition at the horizon location. Depending on the sign of v'{x) we have two different situations. (I) The fluid 
velocity grows monotonically passing from the subsonic to the supersonic region. In this case the external potential 4' 
has a local maximum at the horizon. (II) The fluid velocity decreases monotonically passing from the subsonic to the 
supersonic region. Now the external potential has a local minimum at the sonic horizon. The qualitative behavior of 
the parameters for both cases (I) and (II) is depicted in Fig. ©. 

Also in the case under consideration we can solve the fluid-dynamical equations for a fluid with equation of state 
given by Eq. Q. We have for n 7^ 1 

V \^ J 1 — n \v J 2 

where a, (3 are integration constants. Solving the equation c{v) — v one easily finds the fluid velocity on the horizon, 
Vh = (0^/3""^)^/*^"+^^. Differently from the case of a flux tube with nonconstant section, in this case the formation 
of the horizon does not imply limitations on the range of variation of v. Using the previous equations one can now 
easily show that {dtl) / dv)\v,^ = and {(Pij) / dv'^)\^^ = -{n + 1) so that {d'^i)/dv^)\y^ > (< 0) for n > -1 (n < -1). 
Using the same arguments used for the flux tube with non-constant section, one concludes that for n > —1 {n < —1) 
the horizon forms in correspondence of a maximum (minimum) of the external potential. The value n = — 1 is also 
here a limiting case with the same behavior as that previously discussed. 

In solving Eqs. ^ in the two cases of constant and non constant cross-section, we have not considered n = 1 in the 
equation of state Q). Solutions ||SJ) and © become singular for n = 1 and this case has to be considered separately. 
For n — I the solutions for a flux tube with variable cross-section read 

_ -2 „2 

p(v)—ae 2^, A{v) — — , c = a. (9) 

av 

For a flux tube of constant section we have instead, 

p{v) = 'tp{v) = ~a in — , c = a. (10) 

V V 2 

In both cases we have for the fluid velocity at the horizon = a. A straightforward calculation shows that, as 
expected, the horizon forms for a minimum of the tube cross-section and for a maximum of the external potential. 

As we have already noted one main physical requirement for the formation of an acoustic horizon is that v'{x) 
remain finite and non vanishing throughout the flux tube. In general to keep v'{x) ^Q,oo one needs a fine-tuning of 
the initial conditions. For instance, if v is too large at the entrance of the Laval nozzle the fiow will become supersonic 
at a point upstream from the throat so that v'{x) — > —00 0- This will result in the generation of a shock wave, which 
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will destroy the stationary flow of the fluid. We will not discuss here how this fine-tuning can be realized. We will 
just point out that the condition v'{x) ^ 0,oo can be obtained as necessary conse quen ce of imposing a dynamical 
analogy between the fluid flow and a gravitational, spherically symmetric black hole |20|, |2l| . 

The Einstein equations for a spherically symmetric black hole can be put in correspondence with the fluid motion 
described by Eqs. and constrained by 

2^-X^ + 2e-^^^X'V. (11) 
dr ar ar 

where X — ^(c^ — u^), Y ~ pc, F ^ \n j^, r is the radial coordinate for the black hole spacetime, is the inverse 
of Newton constant and y is a function of the coordinate r, which depends on the particular gravitational model 
under consideration [2ll |. The constraint (|11|) represent a constraint on the geometrical form of the flux tube, which 
forces the cross-section to have a local extremum exactly at the position of the horizon. The constraint Hll|) also 
implies that v'{x) is finite and nonvanishing at the horizon. This can be demonstrated using the trivial identity 
dv/dx = {dv / dY){dY / dX){dX / dr){dr / dx) and showing that each factor in the product is individually finite and 
nonvanishing at the horizon. The third factor is proportional to the Hawking temperature of the horizon, the fourth 
is the fluid density (20L l2l|. therefore both must be ^ 0, oo. The first and the second can be shown to be finite and 
non- vanishing at the horizon by using the explicit solution of the constrained fluid dynamics derived in Ref. |2]| . 

Spherically symmetric flow 

Let us now consider the inviscid, isentropic, spherically symmetric flow of a fluid with equation of state given by 
For simplicity we will set to zero the external forces acting on the fluid. Using spherical coordinates (r, 0, 0), the 
continuity and Euler equations yield, 

^^^=0, .^ + aV-4 = 0. (12) 
dr dr dr 

These equations can be easily generalized to the case when a pointlike source is located at the origin of the spherical 
coordinate system. Whereas the Euler equation remains unchanged the continuity equation gets a contribution 
proportional to S{r). In this more general case for r > the dynamics of the system is still described by Eqs (|12|l . 

The formation of sonic horizons for the spherically symmetric flow can be easily discussed noticing that Eqs. H12|) 
are formally identical to those describing an axi-symmetric flow along a tube with cross-section Air) = Airr"^ . Because 
the section increases monotonically, in view of the previous derived results we should conclude that an horizon cannot 
form. However, these results rely on the assumption that u'(r) = dv/dr is everywhere finite. Relaxing this assumption 
could allow for the formation of an acoustic horizon even though A{r) does not have local extrema. Divergence of v'{r) 
implies that the fluid has inflnite acceleration. The discussion has therefore no direct physical relevance. However, 
considering this more general situation sheds light on the physical mechanism on which horizon formation is based. 

Setting /3 = r^vp, Eqs. (|12|l are solved by Eqs. Q with A{r) = Airr'^. The condition for horizon formation is also 
given by Eq. Using Eqs. (0) and the continuity equation in ifT^ one gets A{v) = 47r^[^^(a — u^)]^/^""^'. One 
can easily show that for n > —1 (n < —1) A{v) has a minimum (maximum) in correspondence of w = w/i, with Vh 
given by Eq. ©. It is immediately evident from the equation dA/dv = {dA/ dr){dr / dv) that in this case local extrema 
of the function A{v) are not related with local extrema of A{r) but with the divergence of dv/dr. The function A{v) 
contains more information about the formation of acoustic horizons then the function A{x) (or A{r)). The vanishing 
of dA/dv is a necessary and sufficient condition for the existence of a point where v — c, . 



* Electronic address: mariano.cadoni@ca. infn.'it1 
^ Electronic address: paolo.pani@ca.infn.iTl 
[1] W. G. Unruh, Phys. Rev. Lett. 46 (1981) 1351. 

[2] H. W. Liepmann and A. Roshko, Elements of Gasdynamics, John Wiley and Sons, 1957. 
[3] M. Visser, arXiv:gr-qc/9311028 

[4] M. Visser, Class. Quant. Grav. 15 (1998) 1767 "arXiv:gr-qc/9712010". 

[5] M. Novello, M. Visser and G. Volovik, "Artificial Black Holes" (World Scientific, 2002) 

[6] M. Visser, C. Barcelo and S. Liberati, Gen. Rel. Grav. 34 (2002) 1719 arXiv:gr-q c7oilllll| . 

[7] V. Cardoso, arXiv:physics/0503042 

[8] C. Barcelo, S. Liberati and M. Visser, arXiv:gr-qc/0505065 

[9] G. E. Volovik, Phys. Rept. 351 (2001) 195 arXiv:gr-qc/000509'T] . 



[10] C. Barcelo, S. Liberati and M. Visser, Int. J. Mod. Phys. A 18 (2003) 3735 arXiv:gr-qc/011003^. 
[11] P. O. Fedichev, and U. R. Fischer, Phys. Rev. A 69 (2004) 033602 arXi v:cond-mat /0303063. . 
[12] G. E. Volovik, Found. Phys. 33 (2003) 349 arXiv:gr-qc/0301043 . ' 
[13] P. O. Fedichev, and U. R. Fischer, Phys. Rev. Lett. 91 (2003) 240407 arXiv:cond-mat/0304342j. 
[14] C. Barcelo, S. Liberati and M. Visser, Int. J. Mod. Phys. D 12 (2003) 1641 arXiv:gr-qc/030506lJ. 
[15] C. Barcelo, S. Liberati, S. Sonego and M. Visser, New J. Phys. 6 (2004) 186 arXiv:gr-qc/0408022]. 
[16] E. Berti, V. Cardoso and J. P. S. Lemos, Phys. Rev. D 70 (2004) 124006 arXiv:gr-qc/0408099 . 
[17] V. Cardoso, J. P. S. Lemos and S. Yoshida, Phys. Rev. D 70 (2004) 124032 arXiv:gr-qc/0410107 . 
[18] H. Nakano, Y. Kurita, K. Ogawa and C. M. Yoo, Phys. Rev. D 71 (2005) 084006 arXiv:gr-qc/041104T1 . 
[19] G. E. Volovik, "The universe in a hehum dropl et", Oxford 2003. 
[20] M. Cadoni, Class. Quant. Grav. 22 (2005) 409 ■arXiv:gr-qc/0410138 



[21] M. Cadoni and S. Mignemi, Phys. Rev. D 72 (2005) 084012 arXiv:gr-qc/0504143l 



[22] S. Liberati, S. Sonego and M. Visser, Class. Quant. Grav. 17 (2000) 2903 arXiv:gr-qc/0003105 
[23] M. A. Sakagami and A. Ohashi, Prog. Theor. Phys. 107 (2002) 1267 arXiv:gr-qc/0108072, . 



